Soft X-ray harmonic comb from relativistic electron spikes 
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We demonstrate a new high-order harmonic generation mechanism reaching the 'water window' 
spectral region in experiments with multi-terawatt femtosecond lasers irradiating gas jets. A few 
hundred harmonic orders are resolved, giving /iJ/sr pulses. Harmonics are collectively emitted by 
an oscillating electron spike formed at the joint of the boundaries of a cavity and bow wave created 
by a relativistically self-focusing laser in underdense plasma. The spike sharpness and stability are 
explained by catastrophe theory. The mechanism is corroborated by particle-in-cell simulations. 

PACS numbers: 52.27.Ny, 42.65.Ky, 52.59.Ye, 52.35.Mw, 42. 65. Re, 52.38.Ph 



High-order harmonic generation, one of the most fun- 
damental effects of nonlinear wave physics, originates 
from many nonlinear it ies, e.g. relativistic [TJ [2] and ion- 
izing matter effects [3]. Harmonics reaching the X-ray 
spectral region enable nanometer and attosecond reso- 
lution in biology, medicine, physics, and their applica- 
tions. Several compact laser-based X-ray sources have 
been implemented, including plasma-based X-ray lasers 
[4], atomic harmonics in gases [3 , nonlinear Thomson 
scattering [5] from plasma electrons [6[ [7] and electron 
beams [8 , betatron radiation [9 , relativistic flying mir- 
rors [10 , and harmonics from solid targets [TJ [2j [TT] 
and electro-optic shocks [12 . Many of them [TJ [5]- 
[TT] come from relativistic laser-matter interactions, de- 
termined by the dimensionless laser pulse amplitude 
ao = eEo/m e coJo > 1 related to the laser peak irradiance 
by Iol = I V e\ol(^/^o) 2 for linear and I oc = 2I 0L for 
circular polarization. Here, e and m e are the electron 
charge and mass, c is the speed of light in vacuum, 
co>o, Ao and Eq are the laser angular frequency, wave- 
length and peak electric field, respectively, and 7 re i = 
1.37 x 10 18 W/cm 2 0. 

In this Letter we experimentally demonstrate bright 
Extreme Ultraviolet (XUV) and soft X-ray harmonics re- 
solved up to a few hundred orders from gas jets irradiated 
by multi-terawatt femtosecond lasers. We propose a new 
mechanism of high-order harmonics generation. 

We have performed two experimental campaigns us- 
ing J-KAREN [13] and Astra Gemini [14J lasers. The 
laser pulse power Po, duration r, wavelength Ao, off-axis 
parabolic mirror f-number, and irradiance in vacuum are 
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FIG. 1. (Color online) The experiment scheme (right) and 
He plasma density profile (left) for the shots in Figs. 2, 3. 



9 TW, 27 fs, 820 nm, //9, and 4 x 10 18 W/cm 2 in the 
first campaign, and 120 TW, 54 fs, 804 nm, //20, and 
4 x 10 18 W/cm 2 in the second campaign, respectively. 
Focusing laser pulses onto a supersonic helium gas jet, 
Fig. [I] we record harmonics in the 80-360 eV spectral re- 
gion in the forward (laser propagation) direction, Figs. [2j 
[3] We use flat-field spectrographs [15] , comprising a 
gold-coated grazing-incidence collecting mirror, spherical 
varied-linespace grating, and back-illuminated Charge 
Coupling Device (CCD), shielded by optical blocking fil- 
ters - two 0.16 fim multilayer Mo/C, two 0.2 fim Pd, or 
one 0.2 /im Ag on 0.1 jam CH substrate. The acceptance 
angle is 3 x 10 -5 sr in the first and 3 x 10 -6 sr in the sec- 
ond campaign. In the latter case, the spectrograph has 
two channels for observation at 0° (forward direction) 
and 0.53° in the laser (linear) polarization plane. The 
wavelength calibration is performed in-place using spec- 
tra of Ar and Ne plasmas. The harmonics energy and 
photon numbers are conservatively estimated using ide- 
alized spectrograph throughputs, i.e. the products of the 
calculated [16] collecting mirror reflectivities and filter 
transmissions, measured diffraction grating efficiencies, 
and manufacturer-provided CCD efficiencies. 
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FIG. 2. (Color) Typical single-shot comb- like spectra, (a) 
Raw data for P = 9 TW, n e = 2.7 x 10 19 cm" 3 , (b) Lineouts 
for two regions marked in (a) by dashed brackets. Dotted ver- 
tical lines for harmonics of the base frequency ujf = 0.885o;o. 
The highest resolved order is n* H — 126. (c) Spectra for the 
shot shown in (a) (red) and for P = 120 TW, n e = 1.9 X 10 19 
cm -3 obtained from the 2 nd (blue) and 3 rd (green) diffraction 
orders. Inset: harmonics in the 2 nd diffraction order; similarly 
for the 3 rd order, (d) The spectral region embracing the 'wa- 
ter window' for P = 9 TW, n e = 4.7 x 10 19 cm" 3 . The drop 
at the carbon K absorption edge, Kc, is due to hydrocarbon 
contamination. Dashed curves in (c,d) for background due to 
the CCD dark current, read-out noise, and scattered photons. 



In both campaigns the laser power significantly ex- 
ceeds the relativistic self-focusing threshold, P s f « 
17 GWx(n cr /n e ), where n cr = racjQ/47re 2 w 1.1 x 
10 21 cm -3 (/im/Ao) 2 is critical density. Thus both lasers' 
pulses self-focus to tighter spots reaching higher irradi- 
ance determined by the plasma density and laser power. 

The comb-like spectra comprising even and odd har- 
monic orders of similar intensity and shape are generated 
by both linearly and circularly polarized pulses in a broad 
range of plasma electron densities from n e ~ 1.7 x 10 19 
to 7 x 10 19 cm -3 . The harmonic base frequency ujf is 
downshifted from the laser frequency c^o, Fig. |2jb,c), in 
correlation with the transmitted laser spectra. 

The data obtained with different lasers demonstrate 
the effect's reproducibility and robustness and the pho- 
ton yield scalability with the laser power, Fig. [2jc) . For 
the 120 TW laser, the photon number and energy in unit 
solid angle within one harmonic at 120 eV are 2 x 10 12 
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FIG. 3. (Color) Modulated spectra with (a) resolved and 
(b) nearly unresolved harmonic structure for Po = 9 TW, 
n e — 4.7 x 10 19 cm -3 . Line-outs (c) and (d) show spectra 
denoted by dashed brackets in (a) and (b), respectively. 



photons/sr and 40/iJ/sr, respectively. These amount to 
4 x 10 9 photons and 90 nJ, assuming 1.5° angular ra- 
dius inferred from our particle-in-cell (PIC) simulations, 
which agrees with similar spectral intensities observed at 
and 0.53° channels. Fig.[2|d) shows the spectrum up to 
the spectrograph throughput cut-off of 360 eV. The emis- 
sion reaches the 'water window' region (284-543 eV), re- 
quired for high-contrast femtosecond bioimaging, where 
the 9 TW laser produces (1. 5+0*4) x 10 10 photons/sr and 
0.8+oJ/iJ/sr, corresponding to 3 x 10 T photons and 1.7 
nJ for the same 1.5° radius. The uncertainties are due 
to the 10% filter thickness tolerance and CCD noise. 

A large number of resolved harmonics, e. g. n* H ~ 
126 in Fig.^b) and > 160 in Fig. [^c, inset), strictly 
bounds the laser frequency change during the harmonic 
emission process, 5uj/uj < l/(2n^), otherwise the orders 
riff and + 1 overlap. The laser frequency decreases 
due to an adiabatic depletion of the laser pulse losing 
energy on plasma waves, while the number of photons 
is conserved [IT7HI9] . For slow depletion, the frequency 
downshift rate equals the energy depletion rate, Suj/uj = 
— 5x/L(±ep [T7tfT9] . For Fig.[2|a,b), the depletion length is 
Ldep ~ 2.7 mm estimated from the observed 70% energy 
transmission through the 0.9 mm plasma. Estimate [T9] 
gives a similar value, -^dep « 8.7(n cr /n e ) 3 / 2 A « 3.4 mm. 
The condition Suj/uj < 0.4% (or, alternatively, the phase 
error of < 25 mrad) gives the harmonic emission length of 
< 12/im. We note that the emission length can be much 
longer than the longitudinal size of the moving source. 

In many shots with linear polarization (40% with the 9 
TW laser), the harmonic spectrum exhibits deep equidis- 
tant modulations, Fig. [3] They are visible with dis- 
cernible high orders and even with nearly unresolved indi- 
vidual harmonics in some shots, Fig. [3^b,d), where blur- 
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FIG. 4. (Color) (a) 3D PIC simulation. Electron density 
n e (top, upper half removed) and its cross-section at z — 
(bottom). The electromagnetic energy density for uj > 4cjo 
(red arcs), (b) Catastrophe theory model. Singularities in 
the electron density created by foldings of the electron phase 
space (x,y, (p y )), where p y is averaged over the laser period. 



ring can be caused by a greater downshift of the laser 
frequency during a longer emission of harmonics, con- 
sistently exhibiting a few times greater photon number. 
Since in some cases the modulation depth is ~100%, we 
conclude that the modulations result from interference 
between two almost identical strongly localized coherent 
sources separated in time and/or space. In shots with 
circularly polarized pulses (90 shots with the 120 TW 
laser), large-scale spectral modulations are not observed. 

The unique properties of the observed harmonics mis- 
match previously suggested mechanisms of high-order 
harmonics generation. Atomic harmonics are excluded 
because both linearly and circularly polarized laser pulses 
produce both even and odd harmonic orders with a weak 
sensitivity to gas pressure. Betatron radiation is not rel- 
evant because the base frequency of its harmonics is de- 
termined by the plasma frequency and electron energy, 
not the laser frequency. Nonlinear Thomson scatter- 
ing cannot provide the observed photon numbers even 
under the most favorable assumptions. For the 9 TW 
shot, Fig. |2ja,b,c), the numerical calculation of single- 
electron radiation spectra [3 [20] gives at 100 eV at most 
2 x 10 -10 nJ/eV sr, assuming the self- focused laser am- 
plitude of ao = 7, pulse duration of 30 fs, and the ob- 
servation angle of 15° close to optimum [7], where the 
high-frequency nonlinear Thomson scattering component 
is most intense. The number of electrons encountered by 
the laser pulse is N e = 7rd^Axn e /4: w 6 x 10 9 , where 
n e = 2.7 x 10 19 cm -3 , d s f « 5/im is the self- focusing chan- 
nel diameter [21], and Ax ~ 12/im is the harmonics emis- 
sion length (see above). This gives at most 1 nJ/eV-sr, 
200 times smaller than experimentally observed. For 
other spectra shown in Fig. [2^c,d) even less fraction cor- 
responds to nonlinear Thomson scattering. 

We performed two- (2D) and three-dimensional (3D) 
PIC simulations of harmonics generation using the code 
REMP [22]. As seen in simulations 231 123], the laser 
pulse undergoes self-focusing [2 , 18 , pushes electrons out 
evacuating a cavity in electron density [I8j [25] and gen- 
erating a bow wave [26]. Fig. [4^ a) represents the case 
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FIG. 5. (Color) 2D PIC simulation, (a) The electron den- 
sity (blue), laser envelope (curves for a = 1,4,7,10), and 
electromagnetic energy density, Wh, for frequencies from 60 
to IOOcjo (red), (b) The upper spike emission spectrum for 
the dashed rectangle in (a), (c) The electron density profile 
10 laser cycles earlier than (a) and the right spike structure 

(d). 



when a linearly polarized (along y axis) laser pulse self- 
focused to the amplitude of ao = 6.6 and FWHM waist 
of 10 Ao, reaches a point where n e = 1.8 x 10 18 cm -3 . The 
simulation window size is 125Ao x 124Ao x 124Ao with 
the resolution of dx = A/32, dy = dz = A/8; the num- 
ber of quasi-particles is 2.3 x 10 10 ; ions are immobile. 
The electromagnetic energy density of a high-frequency 
(uj > 4co>o) field reveals the source of high-order harmon- 
ics, Fig.[4ja, red arcs): they are emitted by the oscillating 
electron density spikes formed at the joining of the bow 
wave and cavity boundaries. While emission is also seen 
from the cavity walls [12] . at higher frequencies much 
stronger radiation originates from the electron spikes. 

The high-order harmonics emitted by the electron 
spikes are seen in 2D PIC simulation, performed in the 
87A x72A window with dx = A /1024, dy = A /112 and 
6 x 10 8 quasi-particles. In Fig. [HJ the laser pulse acquiring 
the amplitude of ao = 10, the duration of 16Ao and the 
waist of 10Ao meets plasma with n e = 1.7 x 10 19 cm -3 . 
Starting from the 7 th order, harmonics are well dis- 
cernible up to the 128 th order, as allowed by the sim- 
ulation resolution. The emission is slightly off- axis with 
the angle decreasing with increasing harmonic order. 

A strong localization of electron spikes, their robust- 
ness to oscillations imposed by the laser and, ultimately, 
their superior contribution to high-order harmonics emis- 
sion is explained by catastrophe theory [27]. The laser 
pulse creates a multi-stream electron flow [28 , 29 stretch- 
ing and folding an initially flat surface formed by elec- 
trons in their phase space, Fig. |4jb), [23] Egg. The surface 
projection onto the (x, y) plane gives the electron density 
where outer and inner folds are mapped into singular 
curves outlining the bow wave and cavity boundaries, re- 
spectively. Catastrophe theory here establishes universal 
structurally stable singularities, insensitive to perturba- 
tions. The bow wave and cavity boundaries produce the 
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'fold' type singularity (A 2 , according to Arnold's classi- 
fication [27]), where the density grows as (A?/) _1//2 with 
decreasing distance to the boundary, Ay. At the joint 
of the two folds, the density grows as (A?/) -2 / 3 , forming 
a higher order singularity - the 'cusp' (A3), Fig. [5|c,d), 
[23| I3T] . Stronger singularities exist [29] , however they 
are not stable against perturbations. Located in a ring 
surrounding the cavity head, the cusp is seen in sim- 
ulations as an electron spike. For linear polarization, 
the harmonic emitting ring breaks up into two segments, 
Fig. Qa). The spike oscillations imposed by the laser 
generate high-order harmonics. The cusp singularity en- 
sures a tight concentration of electric charge, making the 
emission coherent, i. e. the intensity is proportional to 
the particle number squared A 7 " 2 , similarly to the coher- 
ent nonlinear Thomson scattering. However, the cusp 
consists of different particles at every moment of time, 
in contrast to a synchronous motion of the same par- 
ticles. We note that for constructive interference it is 
sufficient that the source size is smaller than the emitted 
wavelength in the direction of observation only. The es- 
timated number of electrons within the singularity ring, 
N e ~ 10 6 , provides a signal level close to the experiment. 

In classical electrodynamics [20 , an oscillating electron 
emits harmonics up to critical order n# c , proportional to 
the cube of the particle energy S e ~ aorac 2 , 

nHc = w c /ujf ~ a*, (1) 

then the spectrum exponentially vanishes as seen in 
our experiments. This gives clues about the harmon- 
ics generation scaling. For the laser dimensionless 
amplitude in the stationary self-focusing [21 , ao, s f = 
(Sir Pon e I P c n cr) 1 / 3 , we obtain the critical harmonic order 
of n Hc ~ Pon e /P c n cr , where P c = 2m 2 c 5 /e 2 « 17 GW. 
The total energy £ s emitted by the electron spike [20] is 
proportional to Af 2 and the harmonic emission time r#: 

£ s « e 2 Ntafocu 2 T H /8c oc N^P^nl^rn . (2) 

Here 7 w (^cr/^e) 1 ^ 2 is the Lorentz factor associated 
with the spike velocity, which is close to the laser pulse 
group velocity in plasma. In the 120 TW shot, Fig. [2|c), 
the energy emitted into a harmonic near 120 eV is 40 
/iJ/sr. In the 2D PIC simulation with the parameters 
close to this shot (Fig. [E]), the energy of 100 th harmonic 
is 30 /iJ/sr. An estimate based on Eq. Q gives 100 
/iJ/sr within the 100 th harmonic, thus providing a simple 
analytical estimate for the expected harmonic energy. 

With the detector at 1.4 m and the source size of 10 jam 
estimated from simulations, the spatial coherence width 
is 1 mm, which is large enough for many phase contrast 
imaging applications in a compact setup [32] . 

In conclusion, irradiating gas jets with multi-terawatt 
lasers we observe comb-like harmonic spectra reaching 
the 'water window' region. We propose the harmonic 
generation mechanism based on self-focusing, nonlinear 



wave formation with electron density spikes, and collec- 
tive radiation by a compact electric charge. Catastrophe 
theory explains the electron spike sharpness and stabil- 
ity. The maximum harmonic order scaling with the laser 
intensity will allow reaching the keV range. Our results 
open the way to a compact coherent X-ray source built on 
a university laboratory scale repetitive laser and accessi- 
ble, replenishable and debris- free gas jet target. This will 
impact many areas requiring a bright X-ray/XUV source 
for pumping, probing, imaging, or attosecond science. 
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